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Abstract 



> 

O ' Nuclear relaxation in half-integer spin chains at low temperatures (T <C J, 

^^ ■ the antiferromagnetic exchange constant) is dominated by dissipation from a 

gas of thermally-excited, overdamped, spinons. The universal low tempera- 
ture dependence of the relaxation rates 1/Ti and I/T2G is computed. 
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Nuclear magnetic resonance experiments have recently been shown to be a powerful tool 
in studying the electronic spin dynamics of the two-dimensional quantum antiferromagnets 
in the cuprate compounds M. The combined measurements of the longitudinal relaxation 
rate, 1/Ti, and the spin-echo decay rate 1/T 2 q, over a wide temperature (T) range allow 
one to learn a great deal about the antiferromagnetic spin- fluctuation spectrum 0-|| . 

Another class of Heisenberg antiferromagnets with novel properties are the one- 
dimensional spin chains with half-integer spins per site. Some nuclear relaxation mea- 
surements on such systems have been performed and more detailed studies are under 
way [0]. Our theoretical understanding of the ground state properties of this system is 
in good shape ||, and there has been limited discussion of the behavior at finite tempera- 
tures ||g|— pTTH . It is the purpose of this paper to use results on the T, wavevector, and frequency 
dependence of the uniform and staggered spin susceptibilities to obtain the NMR relaxation 
rates. There have been earlier discussions [§,13] of the T\ relaxation rate in spin-1/2 chains 
and we will comment on their relationship to our results below. 

An important difference between the two-dimensional antiferromagnets and the half- 
integer spin chains is that the latter are generically critical. This means that the zero- 
temperature spin correlators have power-law decays in both space and time, over a finite 
range of ratios of short-range exchange interactions - this happens because the critical 
fixed-point has no relevant perturbations which respect the symmetry of the underlying 
Hamiltonian. However, for sufficiently large second- neighbor coupling (for example), there 
is a transition to a gapped, dimer phase which is not critical. We will restrict our at- 
tention here to the range of couplings where the ground state is critical. This immedi- 
ately has important consequences for the finite temperature spin-correlators: the entire 
low-temperature region (T <C J) may be considered as the analog of the quantum- critical 
region of two-dimensional antiferromagnets P^||T^ J3[]. There is now no requirement that the 



temperature be larger than some stiffness-associated energy scale, as there is no analog of the 
renormalized-classical []I3] region. Of course, when the small interchain coupling is taken into 



account, three-dimensional long-range order, and the corresponding renormalized-classical 
region, can appear at very low temperatures. 

The principles of conformal invariance allow one to obtain the exact scaling functions 
of the quantum-critical region of many one-dimensional quantum systems [|T5|,|9| P~T]| . Half- 
integer spin chains however posses a complicating feature. There is a marginally irrelevant 
operator at the critical fixed-point, which spoils conformal invariance by inducing logarithmic 
corrections to the leading scaling behavior ]rB|Ji7| . There is no analog of this effect in 
the two-dimensional antiferromagnets. To keep the analysis simple, I will first discuss the 
computation in which this marginal operator is ignored. The logarithmic corrections induced 
by it will be discussed in the next section. 

Let us first look at the staggered spin correlations of the half-integer spin chain. It is 



known fl6| , [TT|| that the equal-time, ground state spin correlators have the staggered compo- 
nent 

(smsAR)) T =o = M-i) fi 0^) 1/2 ^ (i) 

at large spatial separation R. The constant D is non-universal and depends upon the choice 
of microscopic exchange couplings. In the absence of the In R term, this result is sufficient 



to specify the staggered susceptibility at finite temperatures []9| — pLlj] . As noted above, we will 



proceed in the remainder of this section by ignoring the In R term - we will put it back in the 



next section. Then, a simple application of the results of Ref [|TTJ gives us the following result 
for Xs(k,u), the wavevector {k) and frequency dependent staggered susceptibility at finite 
T (for Xs, k is the deviation of the wavevector from ft/a, where a is the nearest-neighbor 
spacing) . 

/I _ . h(u + ckY 

x.s(k, W ) = — 7-n gg+g p <•> *'•• - ~ (2) 

4nk B T 

The quantity c is the T = spinon velocity. Note that the only microscopic input into 
this result for the susceptibility are the values of c and D. A plot of the spectral function 






ImXs(k, lu)/lu was presented in Ref [3 1] for a different value of the critical exponents; here we 
present the imaginary part of the result (^) in Fig [TJ. For small wavevectors, with Tick smaller 
than or around k B T, the peak in the spectral function is at uj = 0, indicating the presence of 
overdamped excitations - spinon excitations interact strongly with other thermally excited 
spinons, acquiring a very short lifetime. We will see below that the NMR relaxation is 
dominated by the contribution of these spinons. At larger k, hck ^> k B T, the peak in the 
spectral function moves to finite uo (See Fig [I]) - these are propagating spinons with a lifetime 
of order Ti/k B T. 

We turn next to the uniform spin susceptibility, Xu(k, to) , where k is now the true wavevec- 
tor, measured from the zone center. Unlike the staggered component, the overall normal- 
ization of the uniform spin susceptibility is not arbitrary, as the total spin is a conserved 
quantity. It is therefore useful to define a magnetization density m a (R) which is the spin 
per unit length. The T = correlator of m is given by || 

(m a (R,T)m b (0,0)} T=Q = --^ , + — — -) (3) 

8n z \{R + icry {R — ict) 1 ) 

where r is the Euclidean time. We can obtain the finite T form of this correlator by 
conformally mapping onto a Matsubara strip, keeping in mind that the magnetization density 
is a component of a current and has non-zero conformal 'spin'. This procedure yields 

/ t-n \ /a nw Sij fk B T\ cosh(27iRk B T/hc)cos(27iTk B T/h) -1 

(mAR, r)m,(0,0)} = — — — ^ (4) 

V V ' n " 2 \ he J (cosh{27rRk B T/hc)-cos{27rTk B T/h)) 2 V ; 

Performing a Fourier transform of this result to wavevectors k and Matsubara frequencies 
u n , we find the remarkably simple result 

Xu ^ lUn) = 2^^FT^I (5) 

Note that, unlike Xs(k,u), all T dependence has disappeared ! The T dependence in (|J) is 
obtained while performing the discrete Fourier transform from frequency to time: the T de- 
pendence is contained entirely in the frequency spacing of the Matsubara sum. Analytically 
continuing to real frequencies we get 



X " ( ^ ) = 2^ c 2 fc 2 -(c, + *e) 2 (6) 

where e is a positive infinitesimal. A conspicuous property of this result is that there is no 
damping of the pole at u = ck, even at finite T. This is, of course, a property only of the 
scaling limit. Upon considering corrections to scaling, some damping should appear, but 
will be suppressed by powers of T/J. Related to the absence of damping, is the fact the 
spectrum of magnetization fluctuations is propagating and not diffusive. The spin diffusion 
constant is effectively infinite in the scaling limit. 

One might, at this point, raise the issue of whether it is legitimate to neglect the damping 
of the uniform magnetization modes: perhaps the damping coefficient will be dangerously 
irrelevant and will contribute a singular T dependence to the relaxation rates computed 
below. This possibility appears to me to be quite unlikely. The only role of temperature in 
all of the computations discussed here is to act as a finite-size cutoff in the imaginary time 
direction to a critical theory: the T = result should surely be obtained in the limit T — ► 0. 

We are now finally in a position to compute the nuclear relaxation rates. We will use 
the following expressions [18J , which are appropriate for the relaxation of a nucleus coupled 



to the electronic spins by a hyperfine term: 



Ti 
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(7) 
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where a is the lattice spacing, and A\\ (A±) are the hyperfine couplings parallel (perpendic- 
ular) to the field; their k dependence is expected to be smooth and arises from appropriate 
form factors. We have also neglected contributions from nucleus-nucleus dipolar couplings 
which could be important in some materials. The susceptibility x should include contribu- 
tions from both the uniform and staggered spin fluctuations. It is now a straightforward 
matter to insert (0) and (|]) into (0) and obtain the T dependence of the rates. Simple 
power counting shows that the contribution of Xs to the rates behaves as 1/Ti ~ T° and 
1/T 2 g ~ T -1 / 2 , while the contribution of Xu scales as 1/Ti ~ T and 1/T 2 g ~ T°. In both 
cases the contribution of the staggered component is dominant for small T. A complete 
calculation of this term yields finally the following leading low T result: 
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T 2 g ~ 2h c \ k B Ta 

where the numerical factor, I is given by the integral 

r((l + zx)/4) 
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71.276591604... (9) 

r((3 + zx)/4) V ; 

We emphasize that these results have neglected the logarithmic corrections to scaling which 
we will consider below. Note that the unknown prefactor D cancels out upon considering the 
ratio of the rates: this should be a convenient way of experimentally testing these results. 



The T independent behavior of 1/Xi (modulo logarithmic corrections discussed below) 
has already been noticed in pervious analyses |||12|]; however a different T-independent value 
was obtained as these works did not account for the damping of the spinon states. 



Logarithmic corrections 

We now consider the consequence of the marginally irrelevant operator present in the 
field theory of half-integer spin chains ||. The basic result is easy to state: both ex- 
pressions for the relaxation rates in @ acquire an identical, multiplicative prefactor of 
In ' (J/T). Further there are subdominant additive corrections which are suppressed by 
powers of l/ln( J/T): the form of these additive corrections will be different for the two 
relaxation rates. As these additive corrections are only logarithmically suppressed, it may 
be necessary to have T significantly smaller than J before the leading results (§) with their 
In 1 ' 2 (J/T) are accurate. 

The arguments for the logarithmic corrections are simple and closely parallel those pre- 
sented in Refs [ IE , I7 |. One begins by writing down the Callan-Symanzik equation for 



l/Ti(T, A), where A ~ J is an ultraviolet cutoff. It is known that this quantity is finite 
in the limit A — > oo after multiplication by a A-dependent renormalization factor Z\. This 
fact can be used to derive a Callan-Symanzik equation for 1/Ti in which the temperature 
T scales under its canonical dimension as it is nothing but an inverse length in Matsubara 
time direction. Integrating the Callan-Symanzik equation (16];|17j to a scale where T ~ A, 
expresses 1/T as In ' (A/T) times the value of 1/Ti in a system in which the coefficient of 
the marginally irrelevant coupling ~ l/ln(A/T). To leading order, we can neglect this cou- 
pling and carry out the latter calculation in the critical theory with no marginal coupling, 
which is exactly what was done above. A similar argument can be made for I/T2G' the 
In ' (A/T) factor will be the same because the rescaling factor Z\ is identical to that for 
1/Ti. However, the perturbative corrections in powers of the coupling constant ~ 1/ ln(A/T) 
should be different in the two rates. 
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FIGURES 

FIG. 1. A plot of the universal spectral weight [(hsT) 2 / D]Imx s {k, io)/hto of the half-integer 
spin-chain as a function of ul = huj/kBT for various values of k = hck/ksT. The values for 
k = 1 have been scaled down by a factor of 1/3. Notice the overdamped peak at k = 1 and the 
propagating peaks at k = 3, 5. 
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